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Abstract 

We show that the effects of neutrino oscillations on 8 B solar neutrinos are 
described well by the first two moments (the average and the variance) of the 
energy distribution of scattered or recoil electrons. For the SuperKamiokande 
and the Sudbury Neutrino Observatory experiments, the differences between 
the moments calculated with oscillations and the standard, no-oscillation mo- 
ments are greater than 3 standard deviations for a significant fraction of the 
neutrino mass-mixing (Am 2 , sin 2 26) parameter space. 

PACS number(s): 26.65.+t, 14.60.Pq, 13.15.+g 
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I. INTRODUCTION 



A full Monte Carlo simulation is necessary in order to understand in detail the results 
of complicated experiments and to estimate their uncertainties. Extensive simulations will 
be especially important for the next generation of solar neutrino experiments JIJ in which 
background effects, energy- dependent sensitivities, and crucial geometrical factors will in- 
fluence the measured rates. Experimentalists will have the time and patience to test a finite 
set of hypotheses against the results of their measurements and massive simulations and will 
report the results to an eagerly- waiting community of physicists and astronomers. 

What additional information will be most useful to report? 

We focus here on the crucial tests of the shape of the continuum neutrino spectrum, 
which are independent of solar physics to a fractional accuracy of (9(10~ 5 ) Qj. Distortions 
of the energy spectra of neutrino-induced electrons could indicate neutrino oscillations. 

We show in this paper that the first two moments of the (recoil or scattering) electron 
energy distribution — the average energy and the variance of the spectrum — can provide a 
compact and informative summary of the expected effects of neutrino oscillations on the 
continuum solar neutrino spectra. It is more practical to report the first two moments of the 
electron energy distribution and their uncertainties than to describe in detail the output of 
an extensive Monte Carlo simulation. Moreover, theorists can conveniently use the measured 
values of the average kinetic energy and its dispersion, instead of a binned energy spectrum, 
to test different models of neutrino interactions. 

We calculate the changes in the first two spectrum moments, that are caused by vacuum 
neutrino oscillations and by matter-enhanced oscillations, in the two-flavor approximation. 
We illustrate our results by realistic numerical studies of the SuperKamiokande H (electron- 
neutrino scattering) and the Sudbury Neutrino Observatory |4]] (SNO, neutrino absorption) 
detectors. 

The structure of our paper is the following. In Sec. II we discuss the basic characteristics 
of the SuperKamiokande and SNO experiments, and present the standard (i.e., no oscilla- 
tion) expectations for the average kinetic energy, (T) , and the variance of the energy distri- 
bution, a 2 . In Sec. Ill we describe our calculation of the spectral moments in the presence 
of neutrino oscillations. In Sec. IV we show the results in the neutrino mass-mixing plane. 
In Sec. V we estimate realistic uncertainties. In Sec. VI we show how SuperKamiokande 
and SNO can discriminate new physics in representative cases. We summarize our work in 
Sec. VII. 

We recall that the SNO experiment will, in addition to the electron spectrum shape 
meaurement, determine the ratio of the charged-to- neutral current interactions of neutrinos 
in deuterium, which is a powerful indicator of neutrino oscillations M. This measurement 
has been investigated in several papers (see the recent works 0-§|) and is not considered 
here. 
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II. BASICS AND STANDARD PREDICTIONS 



In this section we discuss the basic characteristics of the SNO and SuperKamiokande 
experiments, and calculate the standard electron energy spectraQ expected in the absence 
of oscillations. We also evaluate the standard values of the first two moments of the elec- 
tron energy spectra for SNO and SuperKamiokande. We discuss how well the two-moment 
approximation can parametrize the deviations from the standard electron spectra. 



A. Experimental characteristics 

The SNO experiment [Q makes use of a 1 kton heavy-water Cherenkov detector || to 
observe the charged current (CC) neutrino-deuterium reaction: 

v e + d — > p + p + e~ . (1) 

The electron kinetic energy, T, is distributed between and E v — Q, where E v is the 
neutrino energy and Q = 1.442 MeV. The differential cross section for this reaction has 
been recently discussed in detail in ||, to which the reader is referred for specific results and 
additional references |TD| . 



The SuperKamiokande experiment f3[ makes use of a 22 kton (fiducial volume) water- 
Cherenkov detector [O to observe neutrino-electron scattering. The basic reaction is: 

u e + e" -> v' e + e~ , (2) 

in which both charged currents (CC) and neutral currents (NC) contribute. The electron 
kinetic energy T is distributed almost uniformly between and E u /(1 + 2m e /E u ). The 



cross-section is known to first order in the radiative corrections [[T^ ,|10f| . 

If neutrino oscillations u e <-> v x to an active neutrino occur, the pure NC reaction u x -e 
also contributes to the SuperKamiokande signal: 

v x + e" -> v' x + e" , (x = fi, r) (3) 

In the following, we assume x = n for definiteness. 

If the SuperKamiokande and SNO detectors work as expected, high signal-to-noise ratios 
will be achieved above a threshold energy T min ~ 5 MeV. The background assessment is one 
of the most difficult aspects of these experiments. In the absence of reliable knowledge of 
the background, we make the optimistic assumption that there is no significant background 
contamination above a nominal 5 MeV threshold for both experiments. We also assume 
that the detection efficiency is constant above threshold. 

The electron kinetic energy (and direction) will be measured, in both experiments, by 
the Cherenkov light. The distribution of the measured kinetic energy, T, around the true 
kinetic energy, T", can be described by an energy resolution function of the form: 



1 The standard electron energy spectrum at SNO has also been discussed in detail in 
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R(T, T') = —-=exp 



(V - T_+ 6f 



(4) 



In Eq. (f|), the bias term 5 accounts for a possible uncertainty in the absolute energy 
calibration, and the energy-dependent one-sigma width Ay/ scales as \[T' due to the photon 
statistics: 



where A 10 is the energy resolution width at 10 MeV. 



The reference parameters T min , 5, and A 10 used in this work for SuperKamiokande [T3 



and SNO JTj] are given in Table |. In SuperKamiokande, the energy resolution width Aio 
is expected to be measured with 1% accuracy by means of a dedicated, high statistics 
calibration with a tunable linear electron accelerator [13| . 



B. Standard expectations 

We have calculated the standard (i.e, no oscillation) electron spectra in SuperKamiokande 
and SNO, using the detector parameters given in Table 0, the best-estimate cross-sections 
for the reactions in Eq. ([!]) and Eq. (^) \T2\, and the standard 8 B neutrino spectrum 



Hi- 

Table [TJ gives the calculated values of the first two spectral moments of the standard 
distributions, the average electron kinetic (T) and the variance a 2 . 

Figure 1 compares the normalized standard spectrum (Area = 1, solid line labeled STD) 
with three representative spectra calculated assuming neutrino oscillations in vacuum JET] or 
oscillations enhanced by the Mikheyev-Smirnov-Wolfenstein (MSW) effect in the sun.0 
The spectra with oscillations refer to the values of the neutrino mass- mixing parameters Am 2 
and sin 2 26 that best-fit the small-angle MSW (SMA), the large-angle MSW (LMA), and 
the vacuum oscillation (VAC) solutions^ to the solar neutrino problem, given the published 
results of the four pioneering solar neutrino experiments |T9|- p^1 . 

The LMA solution does not induce any appreciable difference from the standard spec- 
trum, since the corresponding electron survival probability, P(E U ) = P{y e —> v e \E u ), is 
almost constant for E u > 5 MeV. Hence the LMA spectrum is almost indistinguishable 
from the standard spectrum in Figs. l(a),(b). However, for the SMA and especially for 
the VAC cases illustrated, there are significant spectral distortions with respect to the no- 
oscillation (STD) case. The deviations in the spectral moments, (T) — (T) and a 2 — a 2 , are 
discussed in Sec. IV. 



2 The MSW regeneration of v e in the earth is not included in the present work. 

3 The best-fit values (Am 2 , sin 2 26) for the SMA, LMA, and VAC cases are respectively: (5.4 x 
1(T 6 eV 2 , 7.9 x to- 3 ), (1.7 x 1(T 5 eV 2 , 0.69), and (6.0 x lO" 11 eV 2 , 0.96) §. 
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C. Two- moment approximation 

In this section we discuss the extent to which the deviations of the first two moments 
characterize the deformations of the electron spectrum induced by neutrino oscillations. 

For the moderately small changes in the standard electron spectrum expected on the 
basis of many oscillation scenarios, the difference between the spectrum with and without 
oscillations (fosc(T) and f(T) respectively) can be approximated by a second order expan- 
sion in T to a good approximation. Defining £(T) as the difference between the spectrum 
after oscillations and the standard spectrum, 

£(T) = /oscCO - f(T) , (6) 

we can write 

eco - £(t) es f(T) + 7 r2 " ( r)g~ g ° 2 ) ' (7) 

where £(T) is the best quadratic approximation to £(T). 

In Eq. (0), j3 and 7 are dimensionless parameters, and (T) and a 2 are the first two 
moments (average and variance) of the standard distribution f(T). The functional form 
adopted in Eq. (|7|) ensures that JdT^(T) = and therefore preserves the normalization. 

Within the 0(T 2 ) approximation, the spectral deformations represent a two-parameter 
family of functions, that can be labeled either by the two parameters (/3, 7) or by the 
deviations and of the first two moments of the new distribution from their standard values, 
(T)o and cr 2 , through the correspondence: 

(T) - (T) ~ J dTf(T) £(T; 0, 7) T , (8) 
a 2 -a 2 ^ |rfT/(T)^(T;/3, 7 )T 2 . (9) 

It is preferable to determine the deviations of the spectral moments instead of fitting the 
(/3, 7) parameters, because the moments are well-defined independently of the functional 
form of the shape deformation and of the order of its expansion in powers of T. Additional 
experimental information (if any) on the third or higher moments could in principle over- 
constrain this determination, but is not essential in the two-parameter approximation of 
Eq. ©. 

The MSW small mixing angle solution (SMA) represents an an important case in which 
two-moment approximation is an excellent approximation. We have checked numerically 
that the deviation £(T) = /sma(^) — f(T) between the the SMA spectrum /sma(^) and 
the STD spectrum f(T) (see Fig. 1) can be parametrized accurately as indicated in Eq. (0). 
In the entire SMA region allowed at 95% CL. |7] by the four pioneering solar neutrino 
experiments, the residuals, divided by the maximum spectral value, are small: 

|£(T) — £(T)| < j 1% (Superkamiokande) , , , 

max f(T) ~ \ 3% (SNO) . v ' 
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The large mixing angle (LMA) MSW solution is well fit by the quadratic approximation 
since, in any event, the LMA spectral distortions are very small (see Fig. 1). 

For the best-fit vacuum oscillation solution (VAC), the residual terms of 0(T 3 ) or higher 
in the shape deviations are as large as 10% (5%) of the maximum spectral value for Su- 
perKamiokande (SNO). The magnitude of the residuals can be made smaller, or larger, by 
varying the VAC solution within the region allowed at 95% C.L. by the available data J7|. 
When the residuals are large, the two-moment approximation may result in a loss of infor- 
mation. Fortunately, precision is not crucial in the cases in which the distortions, £(T), are 
large. In those cases, neutrino oscillations provide a distinctive signature even with an error 
of the order of 10% in the large distortions. 

The spectral measurements in the SuperKamiokande and SNO experiment may initially 
involve relatively large uncertainties. These errors will decrease as more events are accumu- 
lated and as the systematic uncertainties are better understood. If the known uncertainties 
are not much smaller than the changes in the spectrum, even a first-order approximation 
[7 ~ in Eq. (0)] can be adequate to analyze the data. The spectral deformations are then 
in one-to-one correspondence either with (3 (see ||23||), or with the deviation of (T) from the 
standard value (T) (deviation which is proportional to /3, see ||). The use of variations of 
(T) to signal spectral deformations was first emphasized in 0(see also ref. [24]). 

In summary, we see that the description of the spectral deformations in terms of two 
parameters (i.e., the deviations of the first two moments) is a useful approach to the analysis 
of the SuperKamiokande and SNO measurements of the electron energy distribution. The 
loss of information is small for the SMA and LMA solutions. For the vacuum oscillation 
cases in which the distortion is not obvious, the first two moments will also be an adequate 
description. 

We shall show below how the neutrino oscillation parameter space can be explored by 
using the deviations of the first two moments from their standard values. 

III. SPECTRAL MOMENTS WITH NEUTRINO OSCILLATIONS 

The calculation of the electron spectra for dense grids of values of the usual neutrino 
mass-mixing parameters, Am 2 and sin 2 29, can be computer-intensive. One has to calculate 
the differential neutrino cross-section, da/dT, smear it with the energy resolution function, 
and integrate the results over the 8 B neutrino spectrum times the v e survival probability P. 
However, if one is interested in just the (first two) spectral moments, only one numerical 
integration over the neutrino energy is needed, with weight functions that are calculated 
once and for all. In this Section we describe this fast method of calculation for the SNO and 
SuperKamiokande experiments. The application to vacuum and matter-enhanced neutrino 
oscillations will be illustrated in Sec. IV. 

A. Moments of the SNO electron spectrum 

Consider a neutrino of energy E v that is absorbed by deuterium, reaction ([I]). The cross 
section &cc{E v ) for producing an electron that has a measured kinetic energy greater than 

-^min IS. 
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acc(^) = [ dT [ dT'R{T, T>) dacc [f; V) . (11) 

where da/dT' is the differential cross-section for the reaction (|I|) [T^J§], and the energy 
resolution function is given in Eq. (|]). The "hat" is intended to remind the reader that <3"cc 
is a smeared and reduced cross section, i.e., the electrons with measured energy below T min 
are not counted. 

The average value of the measured electron kinetic energy above threshold is: 

t{E v ) = a£(E„) f dTT f dT'R(T, T') ^c(^, T') (12) 
Analogously, the average value of the electron kinetic energy squared is: 

f \E V ) = &&(E V ) [ dTT 2 [ dT'R(T, T) ^c(^, T') ^ 

So far we have considered a neutrino of fixed energy E v . If we average the values of 
T(E V ) and T 2 {E V ) over the 8 B neutrino spectrum, \(E U ), and the v e survival probability, 
P(E V \ we obtain: 

/ dE v X(E U ) P{E V ) a cc (E u ) f{E v ) 

(T) = - , (14) 

C dE u \(E u )P(E u )a C c(E u ) 



[ dE v \(E U ) P{E„) a cc (E u ) f 2 (E u ) 
jdE u \(E u )P(E u )o cc (E u ) 

The variance of the measured electron spectrum is defined as: 

a 2 = (T 2 ) - (T) 2 (16) 
The subscript "0" denotes "standard" (i.e., no oscillation) quantities: 

(T)o = (T) p=i , o* =o*^ (17) 



p=i 



B. Moments of the SuperKamiokande electron spectrum 

The calculations of the spectral moments in SuperKamiokande are slightly more compli- 
cated than in SNO, since the neutrino arriving at the detector can be either a u e or a u^, 
interacting with cross sections a e and respectively. The probabilities of these two occur- 
rences are P and 1 — P respectively, where P is the electron neutrino survival probability. 

For a neutrino of energy E u , we need to calculate the average cross section and the 
average electron kinetic energy for a = e, \i: 
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a a {E u ) = j dT j dT'R(T, T') da ^j^ VI , (i 8) 



f a {E v ) = a-\E v ) j dTT j dT 1 R(T, T') da ^^ . (19) 

The expression for T 2 {E V ) is analogous to Eq. ([H]) with the substitution T a (E u ) — > T 2 (E U ). 
Averaging T over the neutrino spectrum, while taking into account the probabilities of 



v e and interactions at the detector, we obtain: 



JdE u XPa e f e + j 


r dE v X(l-P)a fl f ft 


J 


f dE v \Pa e + 


fdE u \(l-P) a M 



(T) = J — r lj . (20) 



For clarity, we have suppressed the ^-dependences in Eq. (|20|). 

The calculation of (T 2 ) is analogous to Eq. (|20|) , with the substitution T a — > T 2 . The 
variance a 2 is defined as in Eq. (|16|). 

IV. RESULTS IN THE MASS-MIXING PLANE (Am, sin 2 26» 2 ) 

From the discussion of the previous section, it follows that the calculation of the spectral 
moments in SNO and SuperKamiokande requires integrating over energy various products 
of the following quantities: \{E V ) ( 8 B neutrino spectrum); o"cc(^), T(E V ), and T 2 (E U ) 
(SNO experiment); a e (E v ), T e (E v ), T 2 (E U ), a^{E v ), T^E,^), a.ndT 2 (E u ) (SuperKamiokande 
experiment). These ingredients can be calculated at representative neutrino energies, given 
the detector parameters (Table |). 

A numerical table of the above quantities can be found in [[TTJ. The reader can find 
there also tables of the last ingredient needed to calculate (T) and a 2 , namely the oscillation 
probability P(E U ), for the best-fit cases SMA, LMA, VAC shown in Fig. 1. Anyone wishing to 
evaluate the sensitivity of SuperKamiokande or SNO to some neutrino model not considered 
in this paper (e.g., resonant magnetic transitions, neutrino decay, violation of the equivalence 
principle, or exotic transitions ||25|| ) can use the data and software in [fH] |. 

We have calculated the values of the spectral moments (T) and a 2 for continuous ranges 
of the neutrino mass-mixing parameters for two-flavor oscillations. In particular, we have 
considered the rectangle (sin 29, Am 2 /eV 2 ) G [1(T 4 , 1] ® [10" 9 , 10~ 3 ], which is relevant to 
MSW oscillations, and the rectangle (sin 2 29, Am 2 /eV 2 ) G [0.4, 1] <g> [3 X 10 -11 , 2 x 10" 10 ], 
which is relevant to vacuum neutrino oscillations. 

Figure 2 shows the calculated results for the SuperKamiokande experiment and Fig. 3 
gives the results for the SNO experiment. 

Figures 2(a),(b) show the fractional differences (%) in the first two moments of the 
SuperKamiokande electron spectrum, ((T) — (T) )/ (T) and (a 2 — <7q)/<7q, for MSW neutrino 
oscillations. Superposed are the regions allowed at 95% C.L. by current fits to solar neutrino 
data [0 (shaded), with thick dots marking the best-fit points. 

The deviations of the moments are larger in two characteristic regions, corresponding 
to the adiabatic MSW branch (horizontal region), and non-adiabatic MSW branch (slanted 
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region). In the adiabatic branch, there is a strong suppression of the high-energy part of 
the 8 B neutrino spectrum. Therefore the mean value of the electron energy also decreases, 
until there are very few electrons above the experimental threshold energy T min (5 MeV 
nominal value). The variance is also reduced by the combination of the depletion of the 
high energies by oscillations and the threshold cut. The fractional spectral variations are so 
large that part of the adiabatic branch is already excluded by the "low statistics" electron 
spectrum measurements at Kamiokande [EB| . For the non-adiabatic branch (relevant for the 
small mixing angle solution) the low-energy part of the 8 B neutrino spectrum is suppressed, 
so that both (T) and a 2 are shifted towards higher values. The changes in the electron 
spectrum are typically smaller [27] than in the adiabatic region. In the large mixing angle 
region, the suppression of the 8 B neutrino spectrum is almost uniform in the energy range of 
interest and there are no significant deviations of the moments from their standard values. 

Figures 2(c),(d) show the fractional differences in the first two moments of the Su- 
perKamiokande electron spectrum for vacuum neutrino oscillations. Also in this case there 
are two regions where the moments are positively or negatively shifted, corresponding to 
the suppression of low-energy and high-energy 8 B neutrinos. In passing from one region to 
another, there are also cases in which the two moment deviations have opposite sign. These 
intermediate cases correspond roughly to a strong suppression of the peak of the 8 B neutrino 
spectrum, and thus a secondary signature should be a low rate of events. The variations 
of the moments can be very strong for large values of the mixing angles. Asymptotically, 
they become zero for Am 2 — > oo (averaged fast oscillations, uniform suppression), and for 
Am 2 — > or sin 2 29 — > (no oscillation). 

Figs. 3(a-c) are analogous to Figs. 2(a-c), but refer to the SNO experiment. In general, 
the deviations of the moments are larger than in SuperKamiokande, since the final state 
electrons from v-d absorption are more closely correlated with the original neutrino energies 
than the electrons from v-e scattering, and thus are more sensitive to oscillations of the 
parent neutrinos. We will recall this point in Sec. VI. 

Figures analogous to Fig. 2(a) and 3(a) were first reported in 0. In |J, the energy 
resolution function was not included, and the electron energy spectrum in v-d absorption 
processes was approximated by a Dirac delta at any given neutrino energy (which is not 
a sufficient approximation, see the discussion in M). Alternative methods for analyzing 



electron spectra have been proposed in [23|,|2"8 



V. ERROR EVALUATION 

In order to test the null hypothesis of no oscillations, we need to compare the theoretical 
spectral moment deviations shown in Figs. 2 and 3 with realistic error estimates for (T)o 
and a 2 . 

Table [IT] lists the standard values of the first two moments, (T) and a 2 . We will now 
estimate the errors affecting these values by varying the input ingredients of our calculations. 
A more realistic assessment of the uncertainties will be possible after the SuperKamiokande 
and SNO experiments have reported on detector performances. 

Table |TTI| shows a preliminary assessment of the error budget for SuperKamiokande and 
SNO. In Table |T|, the errors due to the energy resolution width and energy scale are obtained 
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by varying the parameters Aio and 5 within the la limits given in Table |. The error due to 
the 8 B neutrino spectrum is obtained by repeating the calculation with the ±3cr 8 B spectra, 
A + and A~ |16[, and dividing the total spread by six. The cross section uncertainty for the 



SNO predictions has been estimated by comparing the values obtained with the Kubodera- 
Nozawa cross-section |L5) (default) and with the Ellis-Bahcall cross-section for the 

u-d CC reaction. For SuperKamiokande, the u-e cross-section is known to first order in the 
radiative corrections [Q. The first order contributions (included by default) alter (T) and 
a 2 by 0.1% and 0.4% with respect to the tree-level cross-section ||3~0| . The contributions due 
to second and higher order corrections are expected to be smaller by a factor of order 1/137 
and can thus be neglected. The statistical errors are calculated for a representative case of 
iV = 5000 collected events, which corresponds roughly to 1-2 years of operation, depending 
on the absolute 8 B neutrino flux and the oscillation scenario. The error of (T) is given by 



a 2 /N . The error of a 2 is given by J [/X4 — (cr 2 ) 2 ]/iV where /14 is the fourth moment of 
the spectrum, /i 4 = ((T — (T)) 4 ). As stated earlier, we have not included uncertainties due 
to the backgrounds; these uncertainties may be important. 

Let us consider the error correlations in Table [TTT[ The correlation between the errors on 
(T)o and <Tq induced by the uncertainties in the energy resolution width Aio (see Table H) 
is equal to one in modulus, since both errors depend on the same parameter (Aio). The 
sign of the correlation is + because any small positive shift in (T) tends to make the low- 
energy part of the electron spectrum wider (above threshold) and thus increases a 2 . For 
similar reasons, the errors induced by the energy scale uncertainty, as well as by the neutrino 
spectrum uncertainties, have correlation +1. The cross-section errors are small (negligible 
for SuperKamiokande) and their correlation can be ignored. The correlation of the statistical 



errors is p = 7x3/ [v o~ 2 y — (o~ 2 ) 2 ] [|3l]], where /X3 is the third moment of the electron energy 
spectrum. 

VI. ISO-SIGMA CONTOURS FOR SUPERKAMIOKANDE AND SNO 

In this Section we illustrate the diagnostic power of SuperKamiokande and SNO to reveal 
possible new physics by measuring electron energy spectra. 

Figure 4 shows contours of equal standard deviations (n-sigma ellipses)^ in the plane 
of the (T) and a 2 deviations that were obtained with the help of Table |TJ. The contours 
are centered around the standard expectations (STD). Also shown are the representative 
best-fit points VAC and SMA. The point LMA is very close to STD and is not shown. 

The cross centered at the SMA best-fit point indicates the solution space allowed at 95% 
C.L. |7|] by the solar neutrino data published so far [see Figs. 2(a), 2(b), 3(a), and 3(b)]. 
The deviations in (T) and a 2 for the SMA solution are confined in a relatively small range. 
For vacuum oscillations, the range of deviations spanned by the whole region allowed at 
95% C.L. by present data [Figs. 2(c), 2(d), 3(c), 3(d)] is much larger and is not indicated 



4 The "number of sigmas" n is defined as n = \fx 2 . The probability content of the error ellipses 
is given by the x 2 distribution for two degrees of freedom . 
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in Fig. 4. The statistical significance of the separation between the SMA and STD points 
in Fig. 4 is dominated by the fractional shift in (T) for both SuperKamiokande and SNO. 
This is not surprising, since the SMA neutrino survival probability increases almost linearly 
with energy for E u > 5 MeV; this increase induces deformations of the electron spectra that 
are nearly linear in T and are well represented by a shift in (T), as discussed in Sec. II (see 
also!). 

The best-fit small mixing angle solution is separated by <^ 3a from the standard solution 
for both SuperKamiokande and SNO (see Fig. 4). The discriminatory power of the two 
experiments appears to be comparable for the SMA solution. The estimated total fractional 
errors of (T) and a 2 in SuperKamiokande are about a factor of two smaller than in SNO 
(see Table |TTT|) . However, the purely CC interaction in SNO [Eq. (|l|)] is a more sensitive 
probe of neutrino oscillations than a linear combination of CC and NC interactions, as 
observed in SuperKamiokande [Eqs. (||]) and ©]•[] Moreover, the final electron energy in the 
u e + d — > p + p + e~ absorption interaction is more closely correlated with the initial neutrino 
energy than in v + e~ — > v' + e~ scattering. The SuperKamiokande detector (according to 
the available information about the likely performances of the detectors) may compensate 
a lower sensitivity to neutrino oscillations with a more precise control of the systematics 
related to the electron energy measurement. 

How do the above results depend upon the energy threshold? We have verified by detailed 
calculations that the statistical significance of the SMA deviations in Figs. 4(a) and 4(b) 
decreases by about 0.6<r per 1 MeV increase in the energy threshold T min . These results 
are valid for both the SNO and the SuperKamiokande detectors and include calculations for 
thresholds of 5, 6, and 7 MeV. 

Although SuperKamiokande and SNO are calculated to be about equally sensitive to the 
SMA solution, different new physics scenarios may be more easily tested by one or the other 
of the two experiments. For instance, the SNO experiment appears to be more sensitive 
than SuperKamiokande to the best-fit vacuum solution (VAC point in Fig. 4). 

When both the SuperKamiokande and the SNO electron spectra have been measured, it 
will be important to estimate the correlation of their errors. In fact, the uncertanties affecting 
the spectra measured in these two different experiments are all independent of each other 
except for the uncertainty in the shape of the input neutrino spectrum [|16|], which is common 
to both experiments. The uncertainties in the neutrino spectrum will induce a correlation 
between the errors in the spectrum shapes measured by SNO and SuperKamiokande. For 
the spectral moments and the detector parameters used in this work, the correlation matrix 
between the total errors of Superkamiokande and SNO is shown in Table |6]. The submatrix 
that links SNO and SuperKamiokande quantities has relatively small entries, p = 0.12-0.17, 
indicating that the "cross-talk" induced by the common neutrino spectrum systematics is 
relatively unimportant. However, this relatively small error correlation may increase or 
decrease for more realistic detector parameters and experimental uncertainties, that will be 
determined after the experiments have been operating for some time. 



5 In practice, the separation of CC and NC events in SNO will be affected by experimental uncer- 
tainties. We have ignored this mis-identification error as well as all the backgrounds. 
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VII. SUMMARY 



In this paper we propose that the electron energy spectra that will be measured in 
the SuperKamiokande and SNO experiments be characterized by their first two moments, 
namely the average electron kinetic energy, (T) , and the variance of the energy distribution, 



a 2 . 



We have shown in Sec. II C that the two-moment approximation is accurate for the MSW 
solutions favored by the four pioneering solar neutrino experiments [|T9|-|2"2|| and is a good 
approximation to the distortions induced by vacuum oscillations except for those cases in 
which the distortion is very large and obvious. 

We have presented in Tables || and |T| the standard values, and estimated uncertainties, 
of (T) and a 2 in the absence of oscillations. We have also presented in Figs. 2 and 3 
the fractional deviations of the two spectral moments induced by neutrino oscillations in 
the whole neutrino mass-mixing parameter space relevant to MSW or vacuum two-flavor 
oscillations. The deviations from the standard values can be greater than 3a, as shown in 
Fig. 4, in a large region of the neutrino parameter space, part of which is currently favored 
by the published results of the four pioneering solar neutrino experiments. The statistical 
significance of the SMA deviations for both SNO and SuperKamiokande decrease by about 
0.6a per MeV for energy thresholds in the range between 5 and 7 MeV|] 
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6 After this paper was completed, H. Sobel suggested we investigate what would be the effect of 
reducing the uncertainty in the energy scale, 5, in the SuperKamiokande experiment from 100 
keV to 50 keV. We find that if 5 is reduced to 50 keV the total uncertainty in both moments is 
reduced by about 15% (5000 events assumed) and the significance of the SMA signal is increased by 
about 0.5a (to about 3.5<r). The uncertainty in the 8 B neutrino spectrum would be the dominant 
systematic uncertainty (excluding background effects) if the energy scale error were as small as 50 
keV. 
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TABLES 

TABLE I. Electron energy threshold and resolution parameters adopted in the analysis. Un- 
certainties are at la. 

Detector parameter Symbol SuperKamiokande SNO 

Kinetic energy threshold ^"min 

5 MeV 5 MeV 

Resolution width at 10 MeV A w 1.6(1 ± 0.01) MeV 1.1(1 ± 0.10) MeV 

Energy scale uncertainty 5 ±100 keV ±100 keV 



TABLE II. Standard predictions for the average kinetic energy, (T), and for the variance, a 2 , 
of the electron spectra at SuperKamiokande and SNO in the absence of oscillations. 

Experiment (T) (MeV) a\ (MeV 2 ) 

SuperKamiokande 7.296 3~42 

SNO 7.658 3.04 



TABLE III. Estimated la uncertainties (A) of the standard moments (T)q and <7q, and their 
correlations (p). Possible background effects are ignored. 





SuperKamiokande 




SNO 






A(T) 


Aa 2 




A(T) 


Aal 




la uncertainties 


(MeV) 


(MeV 2 ) 


P 


(MeV) 


(MeV 2 ) 


P 


Energy resolution width 


0.004 


0.02 


±1 


0.024 


0.10 


±1 


Energy scale 


0.025 


0.06 


±1 


0.052 


0.08 


±1 


8 B neutrino spectrum 


0.016 


0.04 


±1 


0.029 


0.05 


±1 


Cross section 


~ 


~ 


~ 


0.011 


0.01 


~ 


Statistics (N = 5000) 


0.026 


0.08 


±0.63 


0.025 


0.06 


±0.45 


Total la error 


0.040 


0.11 


±0.80 


0.070 


0.15 


±0.83 


Total fractional error 


0.54% 


3.2% 


±0.80 


0.91% 


4.9% 


±0.83 



TABLE IV. The correlation matrix of the total errors affecting the measurements of (T)o and 
Oq at SuperKamiokande (SK) and SNO. The correlations between the SK and SNO entries arise 
from the common uncertainty in the 8 B neutrino spectral shape. 





<T)o (SK) 


al (SK) 


(T) (SNO) 


al (SNO) 


(T) (SK) 
a 2 (SK) 
(T)o (SNO) 
al (SNO) 


1.00 


0.80 
1.00 


0.17 
0.15 
1.00 


0.13 
0.12 
0.83 
1.00 
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FIGURES 

FIG. 1. Standard spectrum of the electron kinetic energy (STD, solid line) for Su- 
perKamiokande (a) and SNO (b). Also shown are three representative spectra that apply if neu- 
trino oscillations occur: small mixing angle MSW (SMA, dot-dashed), large mixing angle MSW 
(LMA, dashed), and vacuum oscillations (VAC, dotted). All spectra are normalized (Area = 1). 
The oscillation cases correspond to best fits [jjj to the results of the four pioneering solar neutrino 
experiments [19 22fl. 



FIG. 2. SuperKamiokande experiment: fractional deviations (%) of the moments of the 
electron energy spectrum caused by two-flavor neutrino oscillations. The results are shown in 
the mass-mixing plane (sin 29, Am 2 ). The shaded regions are favored by current solar neutrino 
experiments at 95% C.L. [7j, with best-fit points marked by a dot. (a) Deviation of (T), MSW 
oscillations, (b) Deviation of a 2 , MSW oscillations, (c) Deviation of (T), vacuum oscillations, (d) 
Deviation of a 2 , vacuum oscillations. 

FIG. 3. SNO experiment: fractional deviations (%) of the moments of the electron energy 
spectrum caused by two-flavor neutrino oscillations. The results are shown in the mass-mixing 
plane (sin 2 29, Am 2 ). The shaded regions are favored by current solar neutrino experiments at 
95% C.L. with best-fit points marked by a dot. (a) Deviation of (T), MSW oscillations, (b) 
Deviation of a 2 , MSW oscillations, (c) Deviation of (T), vacuum oscillations, (d) Deviation of a 2 , 
vacuum oscillations. 

FIG. 4. Iso-sigma contours in the plane of the fractional deviations of the first two spectral 
moments. Labels are as in Fig. 1. (a) SuperKamiokande experiment, (b) SNO experiment. The 
SMA solution can be distiguished at > 3a from the standard (STD) case by both experiments. 
The crosses allow for variations of the SMA solution within the region favored at 95% C.L. by the 
current experiments. See the text for details. 
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